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Abstract 

The Exact-Exchange (EXX) Kohn-Sham formalism, which treats exchange interactions exactly 
within density-functional theory, is applied to one-dimensional periodic systems. The underlying 
implementation does not rely on specific symmetries of the considered system and can be applied 
to any kind of periodic structure in one to three dimensions. As a test system, irans-polyacetylene, 
both in form of an isolated chain and in the bulk geometry has been investigated. Within the EXX 
scheme, bandstructures and independent particle response functions are calculated and compared 
to experimental data as well as to data calculated by several other methods. Compared to re- 
sults from the local-density approximation, the EXX method leads to an increased value for the 
band gap, in line with similar observations for three-dimensional semiconductors. An inclusion 
of correlation potentials within the local density approximation or generalized gradient approxi- 
mations leads to only negligible effects in the bandstructure. The EXX band gaps are in good 
agreement with experimental data for bulk trans-polyacetylene. Packing effects of the chains in 
bulk ircms-polyacetylene are found to lower the band gap by about 0.5 eV. 
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I. INTRODUCTION 



The bandstructure is the key quantity for studying the electronic structure of solids and 
of periodic systems in general. Presently, in most cases, Kohn-Sham (KS) bandstructures, 
calculated either within the local density approximation (LDA) 1,2 or within some generalized 
gradient approximation (GGA) 1 *^, represent the starting point for such investigations. In 
particular, the GW metho d 4 ' 5 ' 6 ' 7 ' 8 , the currently most popular approach for the calculation 
of quasi-particle bandstructures and band gaps, is usually based on LDA or GGA bandstruc- 
tures. While it is possible to perform fully self-consistent GW calculations, this is usually 
not done, because such calculations are not only computationally more demanding than the 
non-self-consistent treatment based on LDA or GGA one-particle states and eigenvalues, but 
also in the few cases where such calculations have been carried out for solids worse results 
compared to experiment^ were obtained. While the non-self-consistent GW approach based 
on KS input is somewhat unsatisfactory from a formal point of view, its success neverthe- 
less underlines the central role of the KS bandstructure. Moreover, also the methods to 
investigate optical properties usually start from the KS bandstructure, as for instance the 
Bethe-Salpeter metho d 8 ' 10 ' 11 ' 12 ' 13 or methods based on time-dependent density-functional 
theor y 14 ' 15 ' 16 ' 17 ' 18 ' 19 or time-dependent current-density-functional theory 20 . 

Despite their central role for the investigation of electronic systems LDA and GGA band- 
structures suffer from severe shortcomings. The most important of these is the substantial 
underestimation of the band gaps of semiconductors and insulator o 21 ' 22 ' 23 ' 24 ' 25 ' 26 . This fact 
was often attributed to the KS formalism itself. Indeed, even the exact KS band gap, i.e., 
the gap which would result from the exact KS exchange-correlaton potential, is not iden- 
tical with the physical band gap 21 ' 22 . The two gaps differ by the derivative discontinuity 
of the exchange-correlation potential at integer electron number o 21 ' 22 . As a matter of fact, 
KS orbitals and eigenvalues originally were considered as auxiliary quantities with no or 
little physical meaning. Recently, however, it was shown that KS eigenvalues and their 
differences are physically meaningful and, e.g., are well-defined zeroth order approximations 
for ionization and excitation energie o 27 ' 28 ' 29 ' 30 ' 31 ' 32 as well as for band gaps 3 ^. In order to 
clarify the meaning of these zeroth order approximations it is instructive to consider how 
excitation energies and band gaps are given in density-functional theory (DFT). Excitation 
energies result from time-dependent DFT. For excitation energies that refer to excitations 
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from an occupied to an unoccupied orbital within a one-particle picture, TDDFT excita- 
tion energies can be considered as the sum of the corresponding eigenvalue difference plus a 
TDDFT correction. Band gaps, on the other hand, are given by the smallest eigenvalue dif- 
ference between occupied and unoccupied orbitals, i.e., the KS band gap, plus the derivative 
discontinuity. At the level of the zeroth order approximation of KS eigenvalue differences, 
TDDFT corrections or derivative discontinuities, respectively, are neglected. Therefore opti- 
cal absorption onsets, i.e., optical band gaps, related to excitation energies and quasiparticle 
band gaps including fundamental band gaps are not distinguished at the level of KS eigen- 
values differences. Nevertheless, eigenvalue difference can be reasonable approximations for 
quasiparticle band gaps and excitation energies as discussed below. 

The shortcomings of LDA and GGA bandstructures, at least in most cases, seem to 
predominantely originate from the exchange-correlation functionals used, rather than from 
limitations of the KS formalism itself. The major problem of the LDA and GGA is the 
incomplete cancellation of the unphysical self-interaction of each electron with itself, which 
is contained in the Coulomb energy and potential. These Coulomb self-interactions raise the 
KS eigenvalues. The more localized valence bands are affected more strongly than the less 
localized conduction bands. As a result of the different impact of the self-interaction error on 
valence and conduction bands the LDA or GGA band gaps are artificially decreased. While, 
in principle, the self-interaction contained in the Coulomb potential is exactly cancelled by 
the KS exchange potential, the LDA and GGA for the latter not nearly lead to a complete 
cancellation. 

Recently, a KS approach for solids treating the exchange energy and potential exactly 
was introduce d 341351361 . This exact-exchange (EXX) approach does not suffer from the 
spurious self-interaction and yields band gaps that are much larger than those from LDA 
or GGA calculations. For standard semiconductors band gaps in excellent agreement with 
experiment were obtaine d 3 ^ 3 'i' i3 ^ 3 ^ i4(|) . In these calculations the correlation potential was 
either completely neglected or approximated within the LDA or GGAs. It turned out that 
the choice of the treatment of the correlation potential has little effect on the bandstructure. 
The agreement between experimental and EXX band gaps for standard semiconductors 
shows that for these materials the approximations made in the EXX approach, the neglect 
of the derivative discontinuity and the approximate treatment of the correlation potential 
together have little effect. For these materials the EXX approach therefore is an alternative 
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to the GW method. The EXX approach, from a formal point of view, has the advantage 
that it is a pure KS procedure, which is carried out self-consistently. Furthermore, the more 
physical EXX band structures can be used instead of LDA or GGA bandstructures as basis 
of G W 41 i 42 , Bethe-Salpeter, or time-dependent density-functional methods^. 

Note that EXX band structures, even in the exchange-only completely different 

from Hartree-Fock (HF) band structures. Both exact exchange-only KS and HF methods 
treat exchange exactly, however, the underlying exchange energy and potential are defined 
differently. The exchange energy in both cases is the exchange energy of the corresponding 
Slater determinant, the HF determinant or the KS determinant, respectively. However, 
the latter two determinants, and therefore also the orbitals building them, are different 
and subsequently also the HF and exchange-only KS exchange energies are different. The 
HF and KS orbitals and thus the corresponding determinants are different because the HF 
and exchange-only KS one-particle equations contain a different exchange potential. The 
KS exchange potential is defined as the functional derivative of the exchange energy with 
respect to the electron density and therefore, by definition, is a local multiplicative potential. 
The local multiplicative KS exchange potential acts in the same way on all orbitals, be they 
occupied or unoccupied. The exact exchange potential therefore cancels the unphysical 
Coulomb self-interactions contained in the Coulomb potential for all orbitals. The nonlocal 
HF exchange potential acts in a different way on occupied and unoccupied orbitals and as 
a result only the occupied orbitals are unaffected by Coulomb self-interactions and describe 
electrons in the system under consideration. The unoccupied orbitals, on the other hand, 
describe electrons in the corresponding iV + 1-particle system with N being the number of 
electrons. As a result HF band structures usually exhibit a much too large band gap and in 
case of metals a singularity at the Fermi level. Exact exchange-only KS bandstructures, on 
the other hand, are not affected by such problems. 

While EXX methods yield very good results for standard semiconductors, an EXX study 
of solid noble gases leads to a somewhat mixed picture 4 ^. The calculated band gaps are 
clearly superior to the LDA values, but are still significantly below the experimental data. 
So far, EXX methods were only applied to three-dimensional periodic systems (as well as 
to atoms 44 and molecule a 45 ^ 46 ). It seems desirable to investigate the performance of EXX 
methods also for other types of periodic systems. In this work the EXX approach is applied 
to the electronic structure of trans-p olyacetylene, as a typical representative of an organic 
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one-dimensional periodic system. One-dimensional organic oligomeres or polymeres are of 
great interest due to their high potential as active materials in new optoelectronic devices. 

Trans-polyacetylene is the prototype example for an organic one-dimensional periodic 
system. While an isolated trcms-polyacetylene chain may be considered as a molecular 
wire or one-dimensional semiconductor, in practice usually fibers of trans-polyacetylene or 
crystalline bulk trans-polyacetylene are encountere d 47 ! 48 . Despite the fact that crystalline 
trans-polyacetylene is a three-dimensional periodic material its properties are mostly deter- 
mined by the periodicity along the polyacetylene chain. We therefore occasionally address 
all forms of trans-polyacetylene as one-dimensional periodic materials in this work despite 
the fact that it is shown below that there exists some nonnegligible packing effects in crys- 
talline trans-polyacetylene. Because trcms-polyacetylene is a prototypical one-dimensional 
periodic organic system it is not surprising that a number of theoretical studies of the op- 
tical properties of this system have been carried Ou t 49 i 50 i 51 i 52 i 53 i 54 i 55 . In Ref. |53j the band 
gap and the optical spectrum of an isolated trans-polyacetylene chain were calculated with 
the GW and the Bethe-Salpeter method, respectively. The results, e.g., the band gap of 2.1 



eV and the absorption spectrum with a single peak at 1.7 eV, agree we 



dat a 5Sl57 , which, however, refer to bulk trcms-polyacetylene. In Ref. 
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1 with experimental 
on the other hand, 



a much larger band gap of 4.1 eV for an isolated trans-polyacetylene chain was obtained 
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by a multireference configuration interaction approach. Finally, in Refs. 
gaps of 3.2 and 3.7 eV, respectively, emerged from M0ller-Plesset perturbation theory for an 
isolated trans-polyacetylene chain. As most of these studies considered isolated chains, but 
relied on experimental data from bulk polyacetylene for comparisons, it seems desirable to 
investigate the influence of the geometry and the arrangement of polyacetylene chains on the 
bandstructure. Such an investigation is the second goal of this work, besides the assessment 
of the performance of the EXX approach for one- dimensional periodic organic systems. 



II. METHOD 

Within the EXX KS approac h 34 i 35 i 36 i 37 i 44 i 58 i 59 i 60 i 61 i 62 , the exchange energy as well as the 
local multiplicative KS exchange potential, not to be confused with the nonlocal Hartree- 
Fock exchange potential, are treated exactly, rather than via some approximate density- 
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functional. For a periodic system the exact exchange energy per unit cell is given by 

E x (r) = -(l/iV k ) £ / /drdr^ :k(r) ^ k '[ r) ^ (r/) ^ k(r/) . (1) 



«/kk' 



/rfr'I s (r,r>,(r') = 2E \W^ L \ck) + c.c.l , (2) 



In Eq. (^), 0t,k are valence orbitals, with crystal momentum k. Orbitals generally are 
assumed to be normalized with respect to the crystal volume. In the considered case of 
non-spin-polarized systems, the occupied orbitals are doubly occupied and spin is taken into 
account by appropriate prefactors of 2. The summation runs over all valence bands, v and 
v', and over all k-points, k and k' within the first Brillouin zone. The number of k-points 
is denoted by iV k 

The KS exchange potential v x is defined as the functional derivative v x (r) = 5E x [p]/5p(r) 
of the exchange energy with respect to the electron density. It obeys the equation 

%k( r )&*(r ) 

which shall be called EXX equation. The ck occuring on the right hand side of the EXX 
equation are unoccupied orbitals forming conduction bands c and the nonlocal operator v^ L 
with the kernel 

^(r,r<) = -£^%^ (3) 

is a non-local exchange operator of the form of the Hartree-Fock exchange operator, but 
built from KS orbitals. The eigenvalues of the orbitals 0„k and ck are denoted by e v ^ and 
e ck , respectively. The independent particle response function X s , i.e., the response function 
of the KS system, is given by 

Xs{ r>) = 5pi*L = 2 E 0:k(r)0 c k(r)0; k (rQ0, k (rQ + ^ 
° v s{*') ~k e^k-e ck 
In Eq. (HJ), v s is the effective KS potential. 

If the exchange potential v x , the response function X s , and the right hand side of the 
EXX equation (j2J) is expanded in an auxiliary basis set, here an auxiliary plane wave basis 
set, then the EXX equation turns into a matrix equation, which can easily be handled 
computationall y 35 ' 3 ^ 37 . 

III. COMPUTATIONAL DETAILS 

Originally the EXX formalism was implemented in a plane-wave code optimized for three- 
dimensional periodic cubic systems 37 by exploiting the cubic symmetry and by using the 
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concept of special k-point o 63 i 64 . For the investigation of non-cubic systems, including one- 
or two-dimensional periodic systems via the super cell ansatz, we programmed a new plane- 
wave implementation, which does not rely on symmetry or special k-points. In super cell 
calculations the directions along the artificial periodicity due to the super cell ansatz are 
treated like the other direction related to real physical periodicity in the sense that no other 
plane wave cutoff's or artificial screenings are introduced. Merely, the number of k-points 
is set to one for those directions in the Brillioun zone that correspond to the directions 
along the artificial periodicity. In order to test that the super cells were chosen large enough 
we additionally carried out calculations with more than one k-point along the directions 
corresponding to the artificial periodicity and checked that results did not change and, in 
particular, that no significant dispersion occured in those directions. 

The cutoff for the plane wave basis set to represent the orbitals and the cutoff for the 
auxiliary plane wave basis set representing the KS response function X a and the exchange 
potential v x in all cases were chosen to be 32 Ry and 12 Ry, respectively. While these cutoffs 
are not sufficient for calculations of total energies, they turned out to be sufficient for the 
calculation of bandstructures. Errors in the resulting band gaps can be estimated to be 
below the order of about 0.05 eV, which is more than sufficient for the purpose of this work. 

Instead of special k-points, we used uniform grids of k-points covering the first Brillouin 
zone. For selfconsistent KS calculations of isolated chains of polyacetylene we used 24 k- 
points in the direction of the chain, which ensured very good convergence with the number of 
k-points. For selfconsistent KS calculations of bulk trans-polyacetylene additional k-points 
were placed in the two directions along the reciprocal unit cell vectors not corresponding to 
the chain direction. Here, a uniform 2x32x2 mesh (32 k-points along the chain direction, 2 
k-points along the two other directions) was employed. The calculations of band structures 
and KS response functions, as usually, are based on the KS potentials obtained in the self- 
consistent calculations, i.e., the band structures and KS response functions were calculated 
by diagonalizing, for the required k-points, KS Hamiltonian operators containing the KS 
potentials from the selfconsistent calculations. (Note that the exact KS exchange potential, 
in contrast to the HF exchange potential, does not depend on the crystal momentum k of the 
orbital it acts upon.) For the KS response functions of isolated chains of polyacetylene we 
used a uniform 1x350x1 mesh of k-points, for bulk trans-polyacetylene a uniform 5x150x5 
mesh was employed. 
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In all calculations all unoccupied states were taken into account for both the construction 
of the response function and the right hand side of the EXX equation. Typical values for 
the number of conduction bands were 3600 (1500) for the isolated chain (bulk) case. 

In the case of EXX calculations the use of pseudopotentials is particularly attractive in 
view of the large number of Fock matrix elements to be evaluated in all-electron calculations. 
For that reason EXX pseudopotentials have been introduced quite earl y 36 i 37 i 65 i 66 . However, 
it has been noticed immediately that a spurious long-range structure develops in the EXX 
pseudopotentia l 36 i 37 i 65 i 66 , if the standard procedure for the construction of normconserving 
pseudopotential o 67l6Sl69 is utilized. This long-range structure, which effectively simulates 
an additional artificial charge on the ion, can be traced to the core-valence interaction^- 
In order to eliminate this structure in a systematic fashion a self-consistent pseudopotential 
construction scheme has been introduced 7 ^, in which a clean asymptotic behavior is enforced 
as constraint. The resulting EXX pseudopotentials have been shown to give very accurate 
atomic excitation energies and spectroscopic constants of molecules. In Ref. it has also 
been demonstrated explicitly that, apart from the long-range structure eliminated via the 
self-consistent pseudopotential construction, the core-valence interaction is less relevant in 
the case of the exact exchange than for the LDA. 

In the present work, all calculations with the exact exchange were based on normconserv- 
ing pseudopotentials obtained by this self-consistent pseudopotential construction scheme. 
In all cases the density-functionals (LDA, exact exchange-only, exact exchange plus approx- 
imate correlation) used for the generation of the pseudopotentials and in the self-consistent 
groundstate calculations of the periodic system were chosen to be identical (pure LDA 
pseudopotentials were generated by the standard procedure^) . For hydrogen the pseudopo- 
tentials contained exclusively an angular momentum component with I = 0. For carbon 
pseudopotentials with s- and p-components were employed, with the s-component chosen as 
local component for all calculations. Additional tests with the p-component as local com- 
ponent for the carbon atom were carried out, but did not lead to significant changes of the 
results. For the cutoff radii values of rf l=0 = 0.9, r^ =0 = 1.2 and r^ l=1 = 1.1 (in units of 
Bohr radii) were used, following Ref.—. 
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IV. RESULTS 



A. Isolated, infinite chain of irans-polyacetylene 

As a first system we consider an isolated, infinite chain of trans-polyacetylene, relying on 
the super cell concept. The polymere is planar and exhibits alternating single and double 
bonds between the carbon atoms, which are accompanied by different bond lengths. The 
unit cell contains a C2H2 unit. The bond distances of the carbon backbone were taken from 
experimental data from X-ray scattering of crystalline trans-polyacetylene^. This geometry, 
which shall be named Gl, has also been used in a recent investigation of electronic correlation 
effects in trans-polyacetylene 55 . The carbon-carbon bond lengths in geometry Gl are given 
by d(C — C) = 1.45A and d(C = C) = 1.36A. The lattice constant along the chain direction 
equals 2.455A4&. A (C — H)-bond length of 1.087A has been found by optimization on the 
HF and MP2 level 50 . The angle Z(C - C, C = C) (with respect to the short (C = C)- double 
bond) has been determined to be 120.057°via HF optimization^ 5 -. 

The system was placed in the the xy-plane with the chain oriented along the ^-direction. 
The unit cell is defined by the vectors a! = (a, 0, 0), a 2 = (0, c, 0) and a 3 = (a/2, 0, 6) with 
the parameter c = 2. 455 A defining the lattice vector along the polyacetylene chain. The 
parameter a was set to a = 8 A. The third lattice vector was chosen in a way that the 
p-orbitals of the carbon atoms of two adjacent chains do not lie on top of each other, but at 
maximal distance. The parameter b was set to b = 9A. The values of a = 8A and b = 9A 
turned out to be large enough that the chains can be considered as isolated. Effects on the 
band gap due to coupling of the chains can be estimated by variation of the super cell size 
to be below 0.05 eV. For further details of the geometry see Appendix. 

The calculated LDA band gap equals 0.78 eV, whereas the EXX calculation yields a band 
gap of 1.63 eV. The bandstructure for these two cases is displayed in Fig. [H adjusted with 
respect to each other by placing the Fermi energy at eV. The LDA and EXX bandstructure 
are qualitatively similar, however, the energetical difference between valence and conduction 
bands is significantly smaller in the case of the LDA, leading to the observed small gap. This 
is in line with the finding that for three-dimensional semiconductors the EXX approach 
increases the commonly too small LDA band gap. Addition of an LDA correlation potential 
to the EXX potential leads to a band gap of 1.65 eV, i.e., increases the gap marginally by 0.01 
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eV, whereas addition of a GGA-correlation potential, here the PBE correlation potential, to 
the EXX potential marginally decreases the gap to 1.61 eV. This shows, again in line with 
the observation for three-dimensional semiconductors, that the treatment of the correlation 
potential has no significant influence on the band gap. Within the EXX scheme, the band 
gap comes close to the value of 2.1 eV resulting from GW calculations 53 , but still stays 
somewhat smaller. The calculated band gaps as well as those from several other methods 
are collected in Table |U 

To investigate the sensitivity of the bandstructure of an isolated polyacetylene chain on 
the geometry, we calculated the EXX bandstructure also in a slightly modified geometry, 



denoted G2, which has also been used in Ref. |53j The geometry G2 is characterized by 
the carbon-carbon bond lengths d(C — C) = 1.44A and d(C = C) = 1.36A 71 . Optimization 
yielded^ 3 a lattice constant along the chain direction of 2.473A. A carbon hydrogen bond 
length d(C = H) = l.lA 53 and an angle Z(C — C, C = C) of 118°have been used (see also 
Appendix). In geometry G2 an EXX band gap of 1.56 eV is obtained, which differs only 
slightly from the result observed for the geometry Gl. 



B. Bulk trans-polyacetylene 

The geometries of the chains building the considered bulk traras-polyacetylene are those 
of the above discussed isolated chains in geometry Gl, i.e., the geometry based on X-ray 
structure^. Within the crystalline t r ans-y olvacet ylen e , each unit cell contains two poly- 
acetylene chains. The unit cell according to Ref. |47| corresponds to a simple monoclinic 
Bravais lattice with the values of a = 4.24A, b = 7.32A and c = 2.46A for the lattice 
constants (see Appendix for further details). 

Again we carried out LDA and EXX calculations, which yielded band gaps at the Ap- 
point of 0.98 eV and 1.64 eV. The corresponding bandstructures are displayed in Fig. Eland 
El However in both cases the band gap at the A-point is not the fundamental band gap. 
Along the line from the T- to the A-point, one observes band gaps (0.81 eV for LDA and 
1.55 eV for EXX), which are reduced in comparison to the gaps at the A-point. However, 
this band gap also is not the fundamental band gap. For LDA as well as EXX calculations, 
the smallest band gaps we found are located at the point (0.5,0.5,0.0), in multiples of the 
reciprocal lattice vectors, (see Appendix for numbering of reciprocal lattice vectors), i.e., at 
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the edge of the Brillouin zone. At this point we calculated band gaps of 0.43 eV for LDA and 
1.18 eV for EXX, which represent, compared to the band gap at the X-point, a remarkable 
dispersion of about 0.55 eV for LDA, resp. 0.46 eV for EXX. This finding is in agreement 
with results from Hartree-Fock calculations which estimated a reduction of the band gap by 
0.6 eV due to interband interactions^. The LDA and EXX bandstructures along the line 
from point (0.5,0.5,0.0) to point (0.5,0.5,0.5) in multiples of the reciprocal lattice vectors are 
shown in Fig. HJ and 

Furthermore, we also checked the influence of the addition of an LDA or PBE correlation 
potential to the EXX potential. Like in the case of an the isolated chain, the correlation 
potential turned out to have an only negligible effect on the bandstructure. 

Optical absorption experiments for trcms-polyacetylene in the considered geometry^ 
give an absorption coefficient rising sharply at 1.4 eV and having a peak at about 1.9 eV. 
A second similar measurement^ gives an identical value for the peak and a value of 1.5 
eV for the onset of the absorption. The lattice constants mentioned within the framework 
of this second measurement^, differ slightly from the geometry of Ref. |5_6j, underlying our 
calculations (a = 4.18A, b = 7.34A, c = 2.455A). Nevertheless, the experimental results 
for the two distinct measurements agree more or less and are in reasonable agreement with 
the band gaps from our EXX calculations. (Remember that no distinction between quasi 
particle and optical band gaps is made at the level of a zeroth order approximation relating 
KS eigenvalue differences to observables.) The data of the band gaps for all discussed cases 
have been compiled in Table ITT1 



C. Reponse functions 

The optical absorption strength obtained by the independent particle, i.e., Kohn-Sham, 
response function for an isolated trans-polyacetylene chain in geometry Gl is shown in Fig. 
IHlfor the EXX-case. The onset of the peak coincides with the calculated fundamental band 
gaps, the absorption maximum lies about 0.1 eV higher than the band gaps. A similar 
behavior was obtained for the Kohn-Sham response function for the LDA-case with the 
corresponding peak shifted to a lower energy. The absorption spectrum resulting from 
the EXX response function is quite similar to the one obtained via the Bethe-Salpeter 
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equation^, which includes electron-hole interactions, i.e., excitonic effects, that are missing 
in an independent particle spectrum. This could be due to a cancellation of two effects. 
For an isolated trans-polyacetylene chain the neglect of electron- hole interactions seems to 
shift — the peak of the absorption spectrum to higher energy values. This could cancel the 
effect that the EXX band gap is smaller than the one obtained within the GW scheme^. 
In any case the similarity between the Bethe-Salpeter and the EXX absorption spectra 
suggests that, beyond a possible shift of the spectrum, electron-hole interactions have little 
effect on the optical spectrum of isolated polyacetylene chains. The absorption spectrum 
of the isolated chain obtained via the Bethe-Salpeter equation agrees quite well with the 
experimental spectrum observed for bulk trans-polyacetylene. However, in the previous 
subsection we showed that packing effects in the bulk reduce the fundamental band gap by 
about 0.5 eV. Therefore it remains to be seen whether Bethe-Salpeter calculations for bulk 
trcms-polyacetylene would yield an absorption spectrum in agreement with experiment or a 
spectrum shifted to lower energy. 

The optical absorption spectrum resulting from the EXX independent particle response 
function for bulk trcms-polyacetylene is displayed in Fig. With a value of about 1.6 eV, 
the maximum of the response function lies close to the calculated values of the two band gaps 
at the X-point and on the path between the X- and the T-point (see Table ITTj) . Compared 
to the Kohn-Sham response function for the case of an isolated, infinite chain of trans- 
polyacetylene, the displayed result for the bulk shows a broader peak. This widened shape, 
as well as the maximum position at about 1.6 eV is in agreement with the experimental 
data for bulk trans-polyacetylene^^. 

V. CONCLUSIONS 

The exact exchange KS approach, implemented within a plane wave pseudopotential 
framework, was shown to yield band gaps for a typical one-dimensional periodic organic poly- 
mere, which agree well with experimental data and with data obtained from GW methods 
and methods based on the Bethe-Salpeter equation. Similar to the case of three-dimensional 
periodic semiconductors, EXX band gaps are substantially larger than their LDA or GGA 
counterparts and thus closer to experiment. The EXX approach therefore promises to be a 
viable tool for the investigation of the electronic structure of organic polymeres, a class of 
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materials with high potential as active compound in new optoelectronic devices. 

The findings of this work suggest that also for polyacetylene EXX band gaps are close 
to experimental band gaps, although in this case no direct measurements of fundamental 
band gaps seem to be available and only comparisons between the EXX band gap and the 
optical band gap of bulk irans-polyacetylene can be made. EXX band gaps differ from the 
exact gaps only by the derivative discontinuity of the exchange-correlation potential and 
the contribution of the correlation potential. As in the case of semiconductors, the inclu- 
sion of an LDA or GGA correlation potential has little impact on the bandstructure, which 
suggests that quite generally the effect of the correlation potential can be considered to be 
negligible. Taken this for granted and assuming that fundamental and optical band gap of 
bulk trcms-polyacetylene are not grossly different, the agreement of EXX and experimental 
band gaps implies that also the derivative discontinuity of the exchange-correlation poten- 
tial is small. However, this point needs further investigation, because it could be that all 
presently available approximate correlation potentials are missing features that affect the 
bandstructure. 

An important outcome of the present investigation is the fact that packing effects in bulk 
trcms-polyacetylene have a substantial quantitative impact on the bandstructure: The band 
gap in bulk trans-polyacetylene is about 0.5 eV smaller than the one obtained for isolated 
chains. In future investigations of organic polymeres such packing effects therefore should 
be taken into account. 
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VI. APPENDIX 

In Tables ITTH II VI and |V] displayed in this Appendix the geometric data for all considered 
systems is compiled. Positions (p,q,r) within the Brillouin zone are given in multiples phi + 
qh 2 + rb 3 of reciprocal lattice vectors b 1; b 2 , and b 3 within the main text. The reciprocal 
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lattice vectors bi, b 2 , and b 3 correspond to the lattice vectors ai, a 2 , and a 3 , respectively. 
The vectors a 2 and b 2 are associated with the direction of the polyacetylene chains. 
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TABLE I: Bandgaps (in eV) for an isolated, infinite chain of trans-polyacetylene calculated with 
various density functionals, the EXX case has been calculated for the two different geometries 
Gl and G2. For comparison, values of Hartree-Fock (HF), M0ller-Plesset theory (MP2), mul- 
tireference configuration interaction (MRCI) and quasiparticle calculations (GW) are given as well. 



HF 

MP2 

MRCI 

GW 

LDA(Gl) 

EXX(Gl) 

EXX(G2) 

EXX-VWN(Gl) 

EXX-PBE(Gl) 



a Ref. 



El 



Ref. 



51 



X 



6.06 a 
3.68 a , 3.96 b 
4.11 c 
2.1 d 
0.78 
1.62 
1.56 
1.65 
1.61 



Ref. 



0, d Ref. 



TABLE II: Bandgaps (in eV) for bulk-irans-polyacetylene. Beside the value for the gap at the 
X-point of the Brillouin zone, a smaller gap between the T- and the X-point {X — > V) has 
been observed. The smallest gap was found at the point (0.5, 0.5, 0.0) (in units of the reciprocal 
lattice vectors) lying on the edge of the Brillouin zone. The experimental values of the band gap 
correspond to the energy, where a sharp rise in the absorption spectrum is observed, i.e., to the 
optical gap. 





Ax 




Ajnin 


LDA 


0.98 


0.81 


0.43 


EXX 


1.64 


1.55 


1.18 


EXX-VWN 


1.65 


1.57 


1.21 


EXX-PBE 


1.61 


1.52 


1.17 


Experiment 






1.4°, 1.5 b 



Ref. 



H, 6 Ref. 
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TABLE III: Geometric data for bulk irans-polyacetylene given in atomic units. Realspace unit 
cell vectors are denoted as ai, coordinates of hydrogen (carbon) atoms are denoted as Hj (Q). 





X 


y 


z 


ai 


o.Ulz44Uo(J 


U. U 


0.0 




n ioi/i/ioi r 7 
-U.lzl44zl f 


4. bo /ooyUo 


0.0 




n n 
U. U 


n n 
U. U 


13.832/994/ 


TJ 

Hi 


i ^^oom no 
-1.76389103 


-1.153571655 


-5.367986055 


H 2 


-2.18160819 


-1.165272865 


1.548413675 


H 3 


2.18160821 


1.165272865 


-1.548413665 


H 4 


1.88533319 


-3.484117395 


5.367986055 


Ci 


-3.44604796 


-1.194232545 


-6.546186695 


c 2 


-0.49945126 


-1.124611975 


0.370213055 


c 3 


0.49945126 


1.124611975 


-0.370213055 


c 4 


3.56749012 


-3.443456505 


6.546186695 



TABLE IV: Geometric data for isolated chain of trcms-polyacetylene (geometry Gl) given in atomic 
units. Realspace unit cell vectors are denoted as ai, coordinates of hydrogen (carbon) atoms are 
denoted as Hj (Ci). 



x y z 

ai 15.11781328 0.0 0.0 

a 2 0.0 4.63927896 0.0 

a 3 7.55890664 0.0 17.00753994 

Hi -2.69959124 1.10770723 0.0 

H 2 2.69959124 -1.10770723 0.0 

Ci -0.64546122 1.11114406 0.0 

C 2 0.64546122 -1.11114406 0.0 



18 



TABLE V: Geometric data for isolated chain of trans-polyacetylene (geometry G2) given in atomic 
units. Realspace unit cell vectors are denoted as a;, coordinates of hydrogen (carbon) atoms are 
denoted as Hj (Q). 





ry< 
Jb 


1 1 

y 


Ai 


ai 


15.11781328 


0.0 


0.0 




0.0 


4.67329404 


0.0 


&3 


7.55890664 


0.0 


17.00753994 


Hi 


-2.69868711 


1.24389067 


0.0 


H 2 


2.69850018 


-1.24389067 


0.0 


Ci 


-0.61981899 


1.21122435 


0.0 


c 2 


0.61981899 


-1.21122435 


0.0 
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FIG. 1: Bandstructure for an isolated, infinite chain of polyacetylene for the LDA (dashed lines) 
and the EXX (solid lines) case 
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FIG. 2: Bandstructure for bulk polyacetylene for the EXX case 
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FIG. 4: Dispersion of the band gap for bulk polyacetylene (EXX case) along the line from point 
(0.5,0.5,0.0) to point (0.5,0.5,0.5) in multiples of the reciprocal lattice vectors 
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FIG. 5: Dispersion of the band gap for bulk polyacetylene (LDA case) along the line from point 
(0.5,0.5,0.0) to point (0.5,0.5,0.5) in multiples of the reciprocal lattice vectors 
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FIG. 6: EXX Kohn-Sham response function for an isolated, infinite chain of polyacetylene 
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FIG. 7: EXX Kohn-Sham response function for bulk polyacetylene 
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